y=c(23,5,7,13,16,18,9,11,12)
x=sample(y,size=5,replace=T)
x
n=5;mi=0;vai=0
m=mean(x)
varmle=(sum((x-m)^2))/n
varmle
vmle=var(x)
vmle
for (i in 1:n)
{mi[i]=(sum(x)-x[i])/(n-1)
vai[i]=(sum((x-mi[i])^2)-(x[i]-mi[i])^2)/(n-1)
}
vaim=mean(vai)
mib=mean(mi)
jackm=n*m-(n-1)*mib
jackv=n*varmle-(n-1)*vaim
jack_d=sqrt(jackv)
mla_d=sqrt(varmle)
jack_d
mla_d